Anomalous Diffusion 
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Recent investigations call attention to the dynamics of anomalous diffusion and its connection 
with basic principles of statistical mechanics. We present here a short review of those ideas and 
their implications. 



1. Introduction. - Diffusion is a fundamental problem 
in statistical physics [H, 0, S 01 and much more simple to 
describe than reaction rates [H, @, 0] , for example. From 
the studies of diffusion, it is possible to obtain a direct 
description of important concepts in physics, such as the 
ergodic hypothesis (EH) 
dissipation theorem (FDT) 




10.1 111 a nd the fluctuation- 
Some of the first 



concepts of statistical mechanics, for instance the Gaus- 
sian distribution of particles, were obtained in the study 
of diffusion, which also led promptly to other concepts as 
the Boltzmann equilibrium distribution. In our study of 
diffusion, we use the Mori formalism [3] which has a well 
defined memory function for non-Markovian systems. In 
this context, the Kubo response function [15| is of special 
importance, since it was formulated in terms of correla- 
tion functions which can be obtained by the use of linear 
response theory and are directly connected with exper- 
iments, such as light or neutron scattering (l6| . In this 
article we display some recent results in this area, and 
we point the reader to selected literature directly related 
to this topic. 

2. Anomalous Diffusion. - Anomalies in the traditional 
Einstein diffusion have been the focus of extensive re- 
search in many disciplines [H, E3, El, 0, (U [H ^ . 
In order to address this problem let us consider the gen- 
eralized Langevin equation (GLE) in the form 



er the gen- 



dP(t) 
dt 



U(t - t')P(t')dt' + F(t), 



(1) 



where P is the particle momentum, n(t) is a kernel, or 
memory function, and F(t) is a random force, which ful- 
fills (F(t)} = and the fluctuation-dissipation theorem 
(FDT): 



(F(t)F(t')) = {P 2 ) eq U(t-t'). 



(2) 



The usual manner to study the diffusive dynamics is to 
investigate the mean square displacement of the particles, 
given by 



lim (x 2 (t)} oc t a . 



(3) 



The exponent a classifies the type of diffusion: for a = 1 , 
we have normal diffusion; for < a < 1, subdiffusion; 
and a > 1, superdiffusion. To obtain average values, 
it is important to define the correlation function for the 



dynamical operator P as 

C P {t) = (P(t)P(O)), 



(4) 



where the brackets () denote an ensemble average. One 
can also define the normalized correlation function 



R(t) = Cp(t)/C P (0), 
which obeys the equation 
dR(t) 



dt 



= - [ n{t- t')R{t')dt'. 

Jo 



Its Laplace transform yields 
R(z 



(5) 



(6) 



(7) 



Recently, an important achievement in the dynamics of 
diffusion has been obtained by Morgado et al. |18| ■ They 
have shown that if 



lim II(z) 



cz 



(8) 



where c is a positive nondimcnsional constant, then 

a = l + v. (9) 

Consequently, the behavior of the memory II(z) for 
small z determines the long range behavior of the dif- 
fusion, Eq. ([3]). This result has found many applica- 
tions: the study of molecular motors [HlHiii anomalous 
diffusio n [l3l . I22I 24, [25|, and the dynamics of dipolar 
chains [26| are a few examples. 

In a subsequent work, Costa et al. [l3| have shown that 
the average value of the momentum is 



(P 2 (t)) 



(P 2 



R 2 (t)[(P 2 (0)) - (P 2 ) eq ] 



(10) 



In this equation the average value will be the equilibrium 
value if 



lim R(t) = 0, 

t — >oo 



(11) 



i. e., irreversibility is a necessary and sufficient condition 
for ergodicity to hold in diffusion. For systems which 
violate the condition given in Eq. (|11[) . there will be no 



2 



ergodicity. Using the final value theorem and Eq. (7]) we 
obtain 

lim R(t) = lim zR(z) = lim l — ^. (12) 

i^oo z— tO z^O 1 + CZ Q 

This relation shows that for most of the diffusive regimes, 
< a < 2, the equilibrium condition, Eq. (jTTj) . holds. 
However, for ballistic diffusion, a = 2, it is not fulfilled 



lim R(t) = lim — !— ^ 0. 

t— >oo z^O 1 + C 



(13) 



This means that ballistic diffusion violates ergodicity 
and the fluctuation dissipation theorem [13]. In other 
words, if the ballistic system is not initially equilibrated, 
then it will never reach equilibrium and the final re- 
sult of any measurement will depend on the initial con- 
ditions. In this situation the EH will not be valid. 
For a > 2, no memory of the initial condition is lost 
(R(t — ► oo) = 1) and the process is not diffusive, being an 
activated process for which the GLE does not work [13j. 
The results of this letter apply to all kinds of diffusion, 
< a < 2, described by a GLE independent of the mem- 
ory range. This gives origin to new studies in ballistic 
diffusion ^SH]- 

For any initial distribution of values -P(O), with 
(F(t)P(0)) = 0, it is possible to obtain the temporal 
evolution of the moments of P, 



(P(t)) = (P(0))R(t), 



(14) 



Again, after an infinite time, we expect that (P{t — » 
oo)) = 0. However, for ballistic diffusion, this is not 
the case, as can be seen from Eq. (fT3")) . The astonishing 
result here is the existence of a residual current such as 
in superfluids. 

3. Nonexponential behavior. - From the above results, 
it is quite clear that the correlation for ballistic diffusion 
will not decay exponentially to equilibrium. Besides the 
ballistic case, any anomalous regime will present nonex- 
ponential decay. Even for normal diffusion, a large num- 
ber of relaxations may be nonexponential [l8j |. There 
are a large number of phenomena where the systems 
do not relax immediately to equilibrium. Those phe- 
nomena, usually associated with non-aging, have nonex- 
ponential relaxation and are most commonly described 
by power laws or stretched exponentials. The study of 
anomalous relaxation has produced quite interesting re- 
sults 0, H H EJ ED, S 0, . 

For a system described by a GLE of the form Eq. fTJ , 
the evolution relies on the noise that drives the particles. 
For a harmonic noise UM 



F(t) 



1 



\J p(uj) cos(uit + 4>(ui))duj, (15) 



where p is the noise spectral density, fcs is the Boltzmann 
constant, and <f> is a set of random phases in the range 



< <t> < 2-7T. A systematic study carried on by Vainstein 
et al. 25[ has shown that the spectral density plays a fun- 
damental role in the description of stochastic processes 
as we shall see. First, the memory function H(t) can be 
easily obtained by the use of the FDT, Eq. @, as 



n(i) 



p(u)) cos(ujt)dLu, 



(16) 



in such a way that the average cancels the random terms 
and obviously the memory is a deterministic even func- 
tion. The Laplace transform of the memory n(z) is an 
odd function in z; therefore, R{z) given in Eq. ([7]) is also 
an odd function in z. This implies that by inverting the 
Laplace transform, R(t) is an even function of t 



R(-t) = R(t). 



(17) 



Second, the condition given by Eq. ([5]) to determine 
the exponent a, Eq ((9]), for a spectral density of the form 



becomes 



lim p(u}) ~ uj 13 , 



(3, /3<1; 
1, (3>1. 



(18) 



(19) 



This shows that a noise spectral density in the form of 
a power law can produce only diffusive motion in the 
region < a < 2. Diffusive motion beyond ballistic is 
not allowed, what can be observed by using the Laplace 
transform. 

To study relaxation we need to know R(t), which can 
be calculated analytically in restricted cases, being ob- 
tained numerically most of the times. A recent study 
shows that exponential decay, power laws or even Mittag- 
Leffler functions are particular cases of a more general 
function [25j ] which approximates the decay. Indeed, con- 
sidering time reversal symmetry }16| and the discussion 
above (see Eq. (fTTjO. the correlation function must be 
even and cannot be any of those forms. We shall expose 
here the conditions under which it is possible to obtain 
an approximately exponential decay, what happens in 
certain circumstances for normal diffusion. In this case, 



7 = lim n(z) 

z— >o 



:P(0). 



(20) 



Consequently, the friction in the usual Langevin equation 
is nothing more than the noise spectral density for the 
lower modes. For an arbitrary memory, the system has a 
rich behavior; even for normal diffusion it is possible to 
show the existence of at least three time ranges [25j|. A 
normal diffusion can be obtained using a spectral density 
of the form \M 



2l 



p{uj) = 



, lo <uj s ; 



(21) 



0, LU > U)„. 
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For a broad band noise spectral density, j/ujg < 1, and 
long times t > 7 -1 it is possible to decouple Eq. ([6]) to ob- 
tain an exponential decay of the form R(t) = exp (— ft), 
which will bring the system to equilibrium, that is, 
R(t -> 00) -> 0. 

Most of the experimental situations where anomalous 
relaxation is present arise in complex, nonlinear or far 
from equilibrium structures in which detailed balance 
does not hold. Good examples can be found in super- 
cooled liquids [27| and in glasses H, 29|. Those systems, 
however, apparently do not have any easy analytical so- 
lution. On the another hand, diffusion can present closed 
solution for the main expectation values, and arises as a 
simple laboratory for the discussion of those properties. 

4. Nonlinear dynamics I: Chain Dynamics - As men- 
tioned before [26j |. chain dynamics is a quite interesting 
subject where the application of those ideas may lead to 
important results. There is a particular problem that 
has attracted our attention: the breaking process of a 
chain of N monomers of mass M subject to a strain S. 
The equation of motion for monomer j, j = 1,2, ..N, 
is 35l.l36j 



M- 



dt 2 

where the forces 
Ga -- 



Gj ~ G j+ i 



„r dX i 

~dt 



+ F(t), (22) 



dU(a + S + y) t 
dy 



(23) 



are derived from the interparticle potential U. In this 
situation, it is simpler to treat the random force F(t) as 
delta correlated, i.e., U(t) = 2j5(t). We then define an 
effective potential and consider the system as a one body 
Kramers problem. However, simulations show that the 
breaking rate is around a hundred times smaller than the 
usual Kramers rate [f|. Several approaches were used in 
order to overcome the problem; one of them was to con- 
sider that the collective motion of the chain generates a 
harmonic noise, which is correlated and has an associated 
memory. Those were to be added to Eq. (f2"2"|) in order to 
consider a non-Markovian analysis 37[. This improves 
the results, but does not solve the problem. This is an 
important issue, because the simulation really describes 
a very important experiment. Polymers used as an ad- 
ditive in a turbulent flow have breaking rates which are 
up to 10 -6 smaller than those computed using the simple 
Kramers theory; the simulation [38] , on the other hand, 
is in perfect agreement with the experiments. 

5. Nonlinear dynamics II: Synchronization - Nonlinear 
dynamics plus noise is an explosive mixture with a large 
number of unexpected results. The study of the evolu- 
tion of maps subject to a common noise shows different 
examples of synchronization [39|, |4(J. Besides that, the 
study of Langevin trajectories shows very nice patterns 
both for systems without memory 4l| and with mem- 



6. Final remarks. - Diffusion, one of the simplest phe- 
nomena in physics, is a starting point for the study of 
simple and complex fluids. Many theorems in statistical 
physics and even proper applications of its formalism, 
such as a "simple linear response theory" , rely on the 
crgodic hypothesis and its validation. Again, diffusion 
shows up as a simple way to address the problem. In 
this work we have discussed various regimes of anoma- 
lous diffusion, which are ergodic in the range of expo- 
nents < a < 2, where a defines the asymptotic be- 
havior of the diffusion, Eq. ([3]). For a = 2, we have the 
special ballistic case, for which ergodicity is not valid, 
as we have seen. In recent years, molecular motors have 
been receiving a lot of attention [2^, [23[ because of their 
large potential for pure and applied science. In this sub- 
ject, a discussion such as the one presented in this paper 
may be very useful. Moreover, there are many situations 
that present violation of the EH, particularly in glassy 
systems [3, [29|, and others were the EH holds. For 
example, dynamical simulations and equilibrium statisti- 
cal mechanics were recently used to treat glass transition 
calculations 34| . The agreement found between the two 
approaches is a strong indication of the validity of the 
EH. This is a quite surprising result for such a complex 
system. Disordered systems are still a large universe in 
which to explore the basic assumptions of statistical me- 
chanics. 
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do Distrito Federal (FAPDF). 
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42J. Synchronization of many different phenomena 
arises continuously in the literature [43| . 
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